Double dispersion of the magnetic resonant mode in cuprates 
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The magnetic excitation spectra in the vicinity of the resonant peak, as observed by inelastic neutron scattering 
in cuprates, are studied within the memory-function approach. It is shown that at intermediate doping the 
superconducting gap induces a double dispersion of the peak, with an anisotropy rotated between the downward 
and upward branch. Similar behavior, but with a spin-wave dispersion at higher energies, is obtained for the 
low-doping case assuming a large pairing pseudogap. 
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The magnetic resonant mode, first observed in the super- 
conducting (SC) phase of YBa 2 Cu 3 6+a; (YBCO) d, has 
been in the last decade the subject of numerous studies, with 
the essential information coming from the inelastic neutron 
scattering (INS) experiments USl- It has been found that 
the peak intensity is highest at the commensurate wavevec- 
tor Q = (it, tt), while its frequency u r shifts with doping. 
More recently, detailed studies of the magnetic response in 
the vicinity of the resonant peak revealed several intriguing 
but quite universal features. While in the SC YBCO the 
stronger component of the resonant mode disperses down- 
wards |4], another branch apparently emerging from the same 
peak shows upward dispersion olaLZIl- Similar features have 
been observed in the underdoped YBCO, whereby the upper 
branch evolves into a spin-wave-like mode at higher energies 
§10]. It is quite remarkable that dispersions for various dop- 
ing show quite consistent anisotropic intensity within the q 
plane Jg, ISA with a rotation angle 45° between the upper 
and lower branch. 

On the theory side there appears to be a consensus that the 
resonant peak can be interpreted as a low-energy collective 
antiferromagnetic (AFM) soft mode, becoming undamped (at 
least underdamped) for T < T c due to the onset of the d x 2_ y 2 
SC gap in the electron-hole excitation spectrum 1 10] . Two 
limits of the same scenario seem to be realized. At optimum 
doping and slightly underdoped YBCO the resonant mode is 
weak 1 2], indicating that the collective mode is weakly bound 
excitonic state within the SC gap 03 EH HE El- On the 
other hand, at low doping the dominant part of the intensity 
of spin fluctuations with q = Q are within the resonant peak, 
so the latter one is closer to an undamped AFM paramagnon 
mode fl5ll . 

The downward dispersion of the resonant mode is within 
the random-phase approximation (RPA) and related theories 
for the dynamical spin susceptibility Xq( w ) fiol fl3l [l4l Ecjll 
a natural consequence of the closing of the d x 2_ y 2 SC gap 
towards the nodal direction of the Fermi surface (FS). The 
RPA seems to capture also some upward component (silent 
band) after the disappearance of the downward branch 1 14]. 

In this work we present results of the memory function ap- 
proach to spin dynamics 1 15], which is capable to capture the 
upward dispersion of the resonant mode. Moreover, it is appli- 



cable also in the low doping regime. At the same time, mem- 
ory function representation offers an appropriate framework 
(broader than RPA) for the general discussion of the INS ex- 
periments. Thus, it will be shown that the explanation of col- 
lective mode properties at low doping implies the existence of 
a large SC-like pseudogap. 

The dynamical s pin susceptibility can be generally ex- 
pressed in the form 11511 



Ll> 2 + LjMq(uj) — UJ 2 
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where the 'spin stiffness' ?7 q = — i([5f. q 5^)1} can be eval- 
uated for within models relevant to cuprates 1 15], while the 
'mode frequency' w q = (^q/Xq) 1 ^ 2 1S related to the static 
susceptibility x q = Xq( w = 0)- The latter is a sensitive 
quantity, so we fix it with the fluctuation-dissipation relation 
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whereby the correlation function C q is better known within 
relevant models and also a more restricted quantity, although 
not directly measured via INS so far. Depending on the damp- 
ing function r q (w) = M q (w), Eq. Q is able to deal with the 
overdamped response in the normal state, with the spin-wave 
dispersion at higher energies (at low doping) as well as with 
the resonant peak in the SC. 

Using the method of equations of motion within the t-J 
model it has been shown that the collective spin fluctuations 
decay into electron-hole excitations 11511 . This leads to the 
lowest-order mode-coupling approximation for the damping 
in the normal state 
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where is the effective spin-fermion coupling fl5Tl and 
Av_(lo) = is the single-particle spectral function. Provided 
the existence of 'hot spots' where the FS crosses the AFM 
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zone boundary (being the case for cuprates at low to interme- 
diate doping) we assume that at low-w quasiparticles with dis- 
persion £k and weight Z k can determine the spectral function 
Ak(u>) = ZkS(ui — £k)- This results in a rather constant r q (cj) 
within the normal state at low -uj and at q ~ Q. Although it 
is derived within the specific prototype model, the form of 
Eq. is quite generic for the damping of the collective mag- 
netic mode in a metallic system, since the lowest-energy decay 
processes naturally involve the electron-hole excitations close 
to the FS. It should be noted that similar ex pres sions appear 
also in theories based on the RPA approach P I ill 1411 . For the 
SC phase at T < T c Eq. Q has to be generalized to include 
the anomalous spectral functions full leading to 

TqM ~ ^kq( W kWk+q ~ Vk«k+q) 2 X 

k 

[/(£k) - f(E k - uj)} S(uj - £ k - £k+ q )] , (4) 

where u£ q = w^Z^Z^q/ ijq, while itk,«k are the usual 
BCS coherence amplitudes and E k = \/e k + A k . 

We are interested in the behavior for q ~ Q. So we 
use for simplicity constant Wkq ~ id. Clearly, the most 
relevant region is the vicinity of the 'hot spots' ko where 
FS crosses the AFM zone boundary. The latter depends on 
the effective quasiparticle band, which we take of the form 
£k = —4?7i^7k — 4?72i'7 k — A i where 7k = (cosfc x + cosfc y )/2, 
7 k = cos k x cos k y . In the following, we assume values 
Vi = m = 0.33, t'/t = -0.33 and t ~ 400 meV, corre- 
sponding roughly to hole-doped cuprates, in particular YBCO 
at intermediate doping. The chemical potential fi is chosen 
so that the volume inside the Fermi surface corresponds to 
Luttinger theorem at particular hole concentration Ch, i.e., 
Vf oc 1 — Ch- r/q in Eq. Q is well known from model 
calculations 11511 and quite restricted in range, and we take 
r/ q = 0.5 t. For the SC gap we assume the d x i_ y i form, 
A q = Ao(cos(7a; — cosq y )/2. Thus we end up with few ad- 
justable parameters at chosen c^: the correlation function C q , 
the effective coupling w and the maximum SC gap Ao. 

Intermediate - optimum doping: Within this regime the col- 
lective mode is heavily overdamped in the normal state. The 
indication for the latter is low intensity of the INS in the rel- 
evant low-energy window. Following the sum rule, Eq. Q, 
this can be made compatible only with modest C q < 1. Fur- 
theron we assume the Lorentzian form C q = Cq/ (1 + q 2 / k 2 ) 
where q = q Q. To be specific, we fix for the presented 
case the 'optimum' doping at c/j = 0.15 with Cq = 1.0 
and k ~ 1.25. The SC gap is roughly known from experi- 
ments and we take Ao = 40 meV. The remaining input is the 
coupling w or equivalently Tq within the normal state. We 
have shown in our analytical derivation within the t-J model 
lll5ll that r q cx t, but it is essentially renormalized due to 
AFM spin correlations, in particular at low doping. We can 
stress that for the appearance of the upper resonant branch 
it is crucial that Tq is not too large, as seems to be inher- 
ent within the RPA fill ll4Tl which otherwise yields within 
the intermediate-doping regime formally quite similar expres- 



sions to our Eqs. For results at intermediate doping we 

use below Tq ~ 0.45i. 




Figure 1: Magnetic fluctuations spectra Xqt^) (arbitrary units) at 
intermediate doping Ch =0.15 for momenta: a) along the x direction 
q = g(l, 0), and b) along the zone diagonal q = q(l, 1). 

The spectra in the vicinity of the resonance x'qi 1 ^ ~ w r) 
for q ~ Q have been partly studied in Ref. 1 15]. Presented 
results show besides a pronounced downward dispersion also 
a weaker upward branch. In Fig. 1 we display Xq( w ) f° r mo " 
menta both along the a>axis, q = q(l, 0), and along the zone 
diagonal q = q(l, 1), while in Fig. 2 we present the planar 
q scans of the intensity Xq( w ) at fixed w, as frequently em- 
ployed in presentation of INS results. 

Following observations can be made on the basis of 
Figs. 1,2: a) both presentations clearly reveal two branches 
emerging from the same coherent resonant mode at u r ~ 
41 meV. Intensity plots of both branches within the q plane are 
square-like around AFM Q (see Fig. 2), however with quite 
pronounced anisotropy. b) For the downward branch the in- 
tensities are strongest along the (1,0) direction. This is con- 
sistent with the faster dispersion along the zone diagonal (1,1) 
(see Fig. la) which reduces intensity relative to the (1, 0) di- 
rection and deforms the constant w-scan into a square-like pat- 
tern, c) The development is more sensitive for u > u> r , still 
the situation with the upward branch is just opposite to the 
downward one. The dispersion is stronger along the (1, 0) di- 
rection and consequently the larger intensity is along the (1,1) 
direction, d) Above the damping threshold uj > 2Ao the up- 
ward branch merges into an incoherent response broad both 
in q as well as in uj. Note, however, that the incoherent part 
still exhausts most of the intensity sum rule, Eq. (|2j, even for 

q = Q 
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Figure 2: Normalized intensity plot of Xq( u ) m me 1 plane at inter- 
mediate doping for selected energies uj below and above the resonant 
peak at uj r ~ 41 meV. 



intensity before reaching this q„. The dispersion of J7 q along 
the (1, 0) direction is substantially weaker, as seen in Fig. la, 
which leads to square-like contours and the (1,0) dominated 
anisotropy in Fig. 2. The upper branch in our analysis ap- 
pears as a exciton-like resonance below next thresholds, i.e., 
w u(q) < ^qi* > 1- The condition for such a resonance is that 
the finite damping r q (w tt ) is not too lar ge. The latter seems 
to be the case within the RPA analysis 1 14], where the up- 
per branch does not emerge from the resonant peak. Finally, 
for w > 2Ao the damping r q (aj) is large and quite constant 
and resonant features disappear, with the remaining strongly 
overdamped AFM paramagnon mode. 




Figure 3: XqX^O (arbitrary units) at low doping ch = 0.1 for q = 
9(1,0). 



Let us give some explanation for the behavior of the collec- 
tive mode as observed in Figs. 1,2. At intermediate (near op- 
timum) doping the normal-state damping is large, Tq > ujq, 
and the collective mode is heavily overdamped in the normal 
state. The sharp resonant peak at Q appears due to the step- 
like vanishing damping Tq(ui < S1q) = within the SC 
phase, where £!q = 2Ak* and k* represents the location of 
the 'hot spot' on the FS. Since the damping cut-off is below 
the characteristic 'mode frequency', Op < cjq, the character 
of the resonant mode is excitonic-like 1 1 dt 1 1 1 II . I.e., it appears 
lower but close to Oq and consequently carries only a small 
part of the whole sum r ule,Eq. © [15].' The dispersion of the 
mode, both the downward 1 14, 16] as well as the upward one, 
is intimately related to the properties of the SC gap A^. As 
noted before fill [Til Il5ll . the damping function r q (w) shows 
for q ^ Q several steps, in contrast to a single step at Q. 
Thresholds are determined by the 'hot spot' condition, i.e., by 
processes of zero-energy electron-hole excitations (in the nor- 
mal state) connecting Fermi surfaces ki?i + 2 = q + K 
where kpi are wavevectors on the FS and K are reciprocal 
lattice vectors. Within the SC phase this leads to steps in the 
damping at f2 q = | Ak F1 1 + | Ak F2 1. Away from q = Q there 
are in general four nontrivial 57 q , i — 1,4 with a possible de- 
generacy for q with a higher symmetry in the Brillouin zone. 

The lowest step at J7 q pushes the downward resonant 
branch as o-v(q) < il q . The latter should for q = q(l, 1) 
approach zero at q n = 2ki?„ where kp„ is the nodal point on 
the FS. It is, however, clear from Fig. lb that the branch loses 



Low-doping: At low doping one expects smaller normal- 
state damping Fq(o;) |15] but at the same time larger Cq, 
which generates spin-wave-like dispersion at larger ui, as ob- 
served in INS 1 8, 9]. Both facts also lead to lowering of 
oc Ch flUl . which at the same time corresponds closer 
to the resonance ui r ~ ujq. Moreover, the resonance peak 
exhausts substantial part of the sum rule. To account for a 
weak anisotropy as well as the spin-wave-like dispersion, we 
use in the calculation for C q the form which incorporates an 
incommensurability <5; = (±<5, 0) and (0, ±5) 

4 1 
C q = A^— (5) 

Y V K 2 + 3|q- Q*l 2 

where Q.; = Q + Si. 

In contrast to the intermediate doping, a direct application 
of the damping, Eq. (0, with a single SC gap seems not to 
be sufficient to describe the observed INS results. First, the 
resonant peak in underdoped YBCO remains broad (not res- 
olution limited) for T < T c , only compatible with a finite 
damping persisting in the SC phase |2l|3|,|9j]. Still, there is 
some signature of a double dispersion |9], although the down- 
ward dispersing mode is much less pronounced. Also, INS 
results reveal a drop of intensity for uj < oj c , where u> c < uj r 
can be interpreted as the spin-gap energy scale 001 showing 
up also in NMR relaxation experiments. 

To account for these observations, we generalize at low 
doping the damping function r q (w), Eq. 0, as follows. The 
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Figure 4: Normalized intensity plot of Xqt^) m the Q plane at low 
doping. The resonant peak is at ui r ~ 33 meV. 



SC gap A q is replaced by a large pseudogap A* with the same 
d x 2_ y 2 -character and resulting r*(w). It is well possible that 
the latter behavior persists also above T > T c , but below the 
pseudogap temperature T < T*, well established for under- 
doped cuprates. However, there is an additional damping, 
T q (w) = r*(cj) + T q (cj), whereby T q (w < Q q ) remains 
finite even within the (pseudo) SC gap. Still, there exists a 
lower scale, possibly the coherent SC (or spin) gap, below 
which r q (cj < uj c ) = 0. Since uj c < uj r the dispersion of 
r q (w > lu c ) does not have any significant effect and we as- 
sume further for simplicity a constant T c . 

In Figs. 3,4 we display results for low T ~ 0, correspond- 
ing to low doping Ch ~ 0.1, where following parameters 
have been adopted: Cq = 1.6, k = 5 ~ 0.3, damping 
T c = 18 meV, T* Q = 60 meV, and gaps A^ = 38.5 meV, 
u> c = 10 meV. While there are similarities with the interme- 
diate doping results in Figs. 1,2, there are also evident dif- 
ferences. The resonant peak is broader, but underdamped, 
due to finite T c < uj t . There is a signature of a downward 
branch with the same anisotropy as for the intermediate dop- 
ing, but this branch is much less pronounced and losing fast 
in intensity. The upward dispersion is stronger and transforms 
at larger u> into the usual isotropic AFM spin-wave dispersion 
with uj oc q. 

Let us discuss the correspondence of presented results with 
INS experiments on cuprates. For the intermediate doping our 
results agree with features seen in optimum or slightly un- 
derdoped YBCO O M III '■ a) pronounced downward disper- 
sion with an enhanced intensity along the (0, 1) direction, b) 



broader and less pronounced upward dispersive branch with 
stronger intensity along the (1,1) direction. Both branches 
with quite similar anisotropy we obtain also for low doping 
in Figs. 3,4, whereby the upper one evolves into spin waves 
at higher uj. This seems to be in agreement with underdoped 
YBCO H 0]. However, it should be pointed out that INS 
reveals in the latter case a much less pronounced downward 
branch, in contrast to optimum doping. 

In conclusion, we have shown that the resonant peak dou- 
ble dispersion and its anisotropy are a notrivial consequence 
of the damping r q (o;) in the SC phase, reflecting the <i-wave 
gap A q structure and related thresholds S7 q . In this sense, the 
INS results on resonant peak serve as a very stringent test for 
the mechanism of the collective mode decay and the structure 
of the SC gap. While our basic assumption lll5ll of the decay 
into electron-hole excitations, consistent with other authors 
JIIIIS. does not offer much freedom of interpretation at in- 
termediate doping, there are several open issues at low doping. 
We get a reasonable explanation of experiments only after as- 
suming two distinct energy scales, i.e., the low spin-gap uj c 
and the large d-wave-like pseudogap A*. Without the latter 
it would be quite hard to account for two branches emerging 
from a single resonant feature at ui r and commensurate Q. We 
should point out that a similar quasi-universal development is 
observed even in non-SC cuprates [18]. Thus, experiments 
on magnetic collective mode dispersion, in particular at low 
doping, can get a new insight on the long-standing question of 
low-energy excitations in pseudogap phase of cuprates. 
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